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Abstract 

Recently it has been shown that the gross structure of the bottomonium spectrum is repro- 
duced reasonably well within the non-relativistic boundstate theory based on perturba- 
tive QCD. In that calculation, however, the fine splittings and the S-P level splittings are 
predicted to be considerably narrower than the corresponding experimental values. We 
investigate the bottomonium spectrum within a specific framework based on perturbative 
QCD, which incorporates all the corrections up to 0{a^mh) and 0{agmh), respectively, 
in the computations of the fine splittings and the S-P splittings. We find that the agree- 
ment with the experimental data for the fine splittings improves drastically due to an 
enhancement of the wave functions close to the origin as compared to the Coulomb wave 
functions. The agreement of the S-P splittings with the experimental data also becomes 
better. We find that natural scales of the fine splittings and the S-P splittings are larger 
than those of the boundstates themselves. On the other hand, the predictions of the level 
spacings between consecutive principal quantum numbers depend rather strongly on the 
scale /i of the operator oc Cyi/(mfer^). The agreement of the whole spectrum with the 
experimental data is much better than the previous predictions when ~ 3-4 GeV for 
as{Mz) = 0.1181. There seems to be a phenomenological preference for some suppression 
mechanism for the above operator. 



1 Introduction 



For a long time most successful theoretical approaches to study the heavy quarkonium spectra 
have been those based on various phenomenological potential models. These model approaches 
have been capable not only of reproducing the charmonium and bottomonium spectra to a 
high accuracy, but also of explaining various other properties of heavy quarkonia such as their 
transition rates and decay rates. Through this success, these phenomenological models have es- 
tablished, essentially, that the heavy quarkonium states can be described well as non-relativistic 
boundstate systems; see e.g. for one of the most recent analyses. On the other hand, the 
problem of the phenomenological approaches is that it is difficult to improve the theoretical 
predictions systematically, and that it is difficult to relate the parameters of the models to the 
fundamental theory. 

Recently there has been new progress in explaining these heavy quarkonium spectra within 
the framework of non-relativistic boundstate theory based on perturbative QCD. It has been 
shown that, by incorporating the cancellation of (9(Aqcd) renormalons contained in the pole 
mass and the static QCD potential, the gross structure of the bottomonium spectrum is re- 
produced reasonably well by the perturbative computation of the spectrum up to 0{l/c^) = 
Ola'^rrib) [0, 0]. Furthermore, it has been shown that the static QCD potential, calculated in a 
series expansion in as up to agrees well with typical phenomenological potentials in the 

region relevant to bottomonium spectroscopy, once the (9(Aqcd) renormalon contained in the 
QCD potential is cancelled against that contained in the 6-quark pole mass |^. Since the 
static QCD potential calculated by lattice simulations is consistent with the phenomenological 
potentials in this region 0, all these calculations are consistent with one another. See also 
P|, which made a direct comparison of lattice results and perturbative predictions of the QCD 
potential in a renormalon-subtracted scheme and found the same consistency. 

The key concept which led to these new results from perturbative QCD, is the follow- 
ing. Conventionally, theoretical calculations of the energy levels of a non-relativistic quark- 
antiquark boundstate closely followed that of a QED boundstate such as positronium: it starts 
from the natural picture that, when an electron and a positron are at rest and far apart from 
each other, they tend to be free particles and the total energy of the system is given by the 
sum of the energies of the two particles (pole masses); as the electron and positron approach 
each other, the energy of the system decreases due to the negative potential energy, so that the 
total energy of the boundstate is given as the sum of the pole masses minus the binding energy. 
When the calculation along the same line was applied to the quark-antiquark system, however, 
the perturbative expansion of the boundstate energy turned out to be poorly convergent, due 
to the contributions from infrared (IR) gluons with wave-lengths of order Aqq^. We can re- 
gard this as reflecting the invalidity of the free quark picture when the quark and antiquark 
are far apart from each other. On the other hand, intuitively we expect that there should 
be a way to calculate the boundstate energy in which the contributions of IR gluons can be 
mostly eliminated. This is because when the boundstate size is sufficiently smaller than Aqqj-,, 
IR gluons cannot resolve the colour charges of the constituent particles, so that they decouple 
from this colour-singlet system. Indeed this idea was theoretically validated in the language 
of renormalons and their cancellation |T0[. As a result, the convergence of the perturbative 
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expansion improved dramatically, extending the predictive power of perturbative QCD beyond 
what could be achieved before. 

Nevertheless, there remain some problems regarding the above prediction of the bottomo- 
nium spectrum from perturbative QCD. Among them, especially interesting is the problem that 
the fine structure in the IPj levels as well as the splittings between the 25* and IPj states are 
predicted to be considerably narrower than the corresponding experimental values. One may 
think that the level of agreement of the theoretical prediction with the experimental data is still 
consistent within errors: according to an estimate based on next-to-leading order renormalons, 
each energy level has a theoretical uncertainty of order Aq(-;j3 ■ (a is the typical size of the 
corresponding boundstate) which may be comparable in size to the 25"-!^, splittings and may 
be much larger than the fine splittings in the IPj levels; compare the error estimate in [Q. One 
should note, however, that if we calculate these level splittings instead of the individual energy 
levels, the 0{Aqqy)) renormalons should get largely cancelled when we take the differences of 
the energy levels. Hence, the theoretical uncertainties of the splittings can be much smaller 
than those of the individual energy levels, and the disagreement between the theoretical pre- 
diction and the experimental data may be quite serious. It is important to clarify whether it is 
possible to resolve these disagreements in the fine splittings and the S-P splittings within the 
context of perturbative QCD, e.g. by including higher-order corrections, or whether we need 
to take into account specific non-perturbative contributions for this purpose. 

In this paper we investigate this problem of the fine splittings and the S-P splittings in the 
bottomonium spectrum within a specific framework based on perturbative QCD. We note that 
these splittings have been successfully reproduced by the phenomenological potential models, 
and that a connection between the static QCD potential and phenomenological potentials has 
been elucidated in [^, ^. In order to take advantage of these results, we develop a frame- 
work which enables detailed comparison of the predictions of the phenomenological models and 
of perturbative QCD. We also incorporate some of the higher-order corrections to the non- 
relativistic Hamiltonian of the quark-ant iquark system which have not been included in the 
analysis 0. 

The basic theoretical ingredients of our analysis are as follows: (1) We take into account 
the cancellation of the (9(Aqcd) leading renormalons by reexpressing the 6-quark pole mass 
in terms of the MS mass. (2) We take a specific scheme for the perturbative expansion such 
that all the corrections up to C(a|mb) and 0(a|m6) are incorporated in the calculation of 
the S-P splittings and the fine splittings, respectively. Furthermore, some of the higher-order 
corrections, which appear to be important for these observables, are incorporated. 

The organization of the paper is as follows. In Sec. 2 we present the framework of our calcu- 
lation. We examine the energy levels and the wave functions of our zeroth-order Hamiltonian 
in Sec. 3. The analysis of the fine splittings is given in Sec. 4 and that of the S-P level splittings 
in Sec. 5. Then, in Sec. 6, we compare the whole structure of the bottomonium spectrum given 
by our prediction, by other theoretical predictions and by the experimental data. Concluding 
remarks are given in Sec. 7. We derive a formula useful for our analysis in the Appendix. 
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2 Framework of Calculation 



2.1 Hamiltonian up to C>(l/c^) 

We first recall the non-relativistic Hamiltonian of a quark and antiquark pair given in the 
series expansion in 1/c up to 0{l/(?) which is determined from perturbative QCD. (See e.g. 
| |TT| , |T^, |TB|, [T^O Considering an application to the bottomonium states, we assume the quark 
(antiquark) to be the 6-quark (6-quark). The Hamiltonian is given by 



H = H^ + U + Wa + W, 



NA- 



(1) 



We choose the zeroth-order part of the Hamiltonian to be 



nib 



Ho = 2mb + — + VQCD[r} 



(2) 



where nib is the pole mass of the 6-quark, and Vqcd('^) denotes the static QCD potential up to 
(9(a|). This choice differs from the usual zeroth-order Hamiltonian of the 1/c-expansion, since 
Ho also includes the 0(a|) = 0{l/c) and 0(a|) = 0{l/c'^) terms of the QCD potential. Other 
operators of Eq. (|ID are treated as perturbations to Hq, all of which are 0(l/c^) in the usual 
order counting in 1/c-expansion. U + Wa constitutes the (9(l/c^) part of the Breit Hamiltonian 
known from QED boundstate theory, where the spin-dependent operator is given by 



U = UlsL-S + Us 



U 



LS 



Us 



2nilr^ 



Uo (2^2 
Uo = 



3m^ 



(3) 
(4) 



and the spin-independent operator is given by 



W 



P 



A 



4ml 



+ 



mf 



2m^r 



-2 1 



(5) 



On the other hand. 



W, 



NA 



2mbr'^ 



("i)\2 



(6) 



represents the operator characteristic to the non-abelian gauge theory. In this paper, unless the 
argument is specified explicitly, a^^'^ denotes the strong coupling constant renormalized at the 
renormalization scale n, defined in the MS scheme with ni active flavours, i.e. a^"'^ = ag^'\iJ,); 
Cp = 4/3 and Ca = 3 are the colour factors; L and S are, respectively, the orbital-angular 
momentum and the total spin of the quark-antiquark pair. For the bb system, ui = 4. 



3 




1 2 3 4 5 6 

Figure 1: Construction of the improved potential. 

2.2 Improved Potential 

In our analysis of the bottomonium spectrum and wave functions, we use an improved "poten- 
tial" i?imp(r) instead of 2mb + Vqcd('") in the zeroth-order Hamiltonian Hq. This Ei^p{r) is 
constructed in the following way: We divide the range of r into three regions by introducing 
ultraviolet (UV) and infrared (IR) distance scales, ruv and tir, see Fig. |1|. 

(i) At intermediate distances, ruv < t < rjR, where the bulk of the bottomonium wave 
functions are located, Ei^p{r) is identified with the total energy of the static bb system 
^tot('^)l«(r) = [2?^6 + VQCB{r)]\^,=f,^{r) computed in §. E^^^{r) is defined in Eq. (5) 
of that paper and depends on the parameters (/i), fi, and mc. these are, respec- 
tively, the strong coupling constant defined in the MS scheme with 3 active flavours, the 
renormalization scale, the 6-quark MS mass renormalized at the 6-quark MS mass scale 
and the same for the c-quark. The 6-quark pole mass is reexpressed in terms of fuf,, 
such that the C(Aqcd) renormalons are cancelled in _E^Q^(r). In identifying Ei^p^r) with 
E^^^-{r), the scale /i is determined as a function of r according to the second prescription 
(yU = fi2{r), Eq. (14)) in since E^l!^^{r) is less /^-dependent and its series expansion 
converges better if we choose a larger value for fi when r is smaller, and if we choose a 
smaller value for /i when r is larger, we consider our choice of fi to give a more accurate 
prediction for E^j^^{r) than choosing some fixed (r-independent) value of /x; see Figs. |],^ 
The c-quark MS mass is taken as nic = 1.243 GeV 0. We will explain how we fix frif, in 
our analysis below. For other details, we follow the convention of Sees. IIA and IIB of [^]. 

(ii) At short distances, r < ruv, we use a renormalization-group improved QCD potential. 
It is obtained by integrating the three-loop renormalization-group improved interquark 
force F{r) = —dVQCT){f')/dr, following the method of Sec. 4 of 0]. There, it was shown 



4 



10 



9.5 



> 

o 



Kl 8.5 

a) 













/X = 1 GeV 







r[GeV-^] ^ 




1 r 2 1, 3 

r [GeV-^] 



Figure 2: 



bb I 
totV 



2m6 + VQCD(r) up to O(a^) for = 1, 2, 3 when IJ, is fixed independently of r: 



a) ;u = 1 GeV, and b) = 2 GeV. The scale /U that provides the best convergence of the perturbative 
series depends on r. 




0.6 0.8 1 1.21.41.51. 

//[GeV] 

Figure 3: Dependence of E^^^{r) on the scale /j, for three different values of r = 1, 2, 3 GeV (solid, 
dashed, dash-dotted, respectively). For large r, the flat region (less //-dependent region) moves 
to smaller scales. 
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that the QCD potential becomes more convergent if we improve the interquark force by 
means of the renormahzation group and integrate it over r, rather than directly improving 
the QCD potential by the renormahzation group.Q The initial value for the renormahza- 
tion group evolution of F(r) and the constant part of Ei^p{r) are determined such that 
Eirap{r) becomes continuous at r = ruv up to the first derivative. 

(iii) At long distances, r > tir, we use a linear potential, Ei^p{r) = Ci r + C2, where Ci 
and C2 are determined such that Ei^p^r) becomes continuous at r = riR up to the first 
derivative. 

The main part of the improved potential Ei^p{r) is that given in the intermediate-distance 
region (i). Indeed this part of the potential dictates the main features of the results of our 
analysis. According to its construction, however, i?imp(r) as defined in the region (i) becomes 
unstable and unreliable at r < and r > 1/Aqcd- This is the reason why we continue it to 

other definitions at short and long distances. The short- and intermediate-distance parts (i) 
and (ii) are determined by perturbative QCD, whereas the long-distance part (iii) is not. We 
will show that the shape of Ei^p^r) in the long-distance region (iii) does not affect the results 
of our analysis significantly. It should be noted that, with our definition, the /i-dependence of 
Eimp{r) has been completely eliminated. The /^-dependences mentioned later in this article are 
those which stem from the other terms of the Hamiltonian (see below). We urge the reader 
to consult ^ for the physics background and detailed features of Ei^p{r) as defined in the 
regions (i) and (ii). In our analysis below, we set ruv = 0.5 GeV~^ and riR = 4.5 GeV~^. 



2.3 Perturbative Expansion 

We solve the non-relativistic Schrodinger equation numerically with the improved potential, 

?2 



(imp) 



h -t/j 

nib 



imp \^ ) 1 



(7) 



and determine the zeroth-order energy level E^^ and wave function oi a quarkonium state. 
We treat U, Wa and Wna as perturbations and calculate the perturbative corrections to the 
energy level by 



AE^ = {i;\{U + WA + W,,A)\i^)- 
We will also examine corrections induced by some of the 0{l/c^) operators: 



SU 



5Uls L-S + 5Us 



a 



(ni) 



TX 



1) 



Ca 



55\ 2 



9 



(8) 



(9) 



(10) 



* Here, the renormahzation-group improvement of Vqcd(?') refers to that using the l^-scheme beta function 
(the second paper of 1 12 ) and is different from the renormahzation-group improvement in pO|-p3) (resummation 
of the next-to- next-to-leading logarithms) . 
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SUs 
5Wna 



„("i)N2 




'11^ 



a 



(ni) 



X 



IT 



f3o 



ai 



where 



P 

ai 
b2 



V 36 

•jE = 0.5771... is the Euler constant; Tr 



log(/ir) +7^;, 
log(mbr) + -fE, 

y ~ 3 

^-hog2^CECA-lcETnni 



36 



(12) 



(13) 
(14) 

(15) 
(16) 



(17) 



1/2. The spin-dependent operators 6Uls and 6Us 

We derived 



(Earher incomplete results can be found in [0, |Tl 



were derived in [|1^, [16 

SWjsfA from the result of [|T^] in the following way: We discard the logarithm originating from 
the IR divergence associated with the energy scale. This gives the first term of Eq. (p!^). 
We have to take into account the contributions which come from the unitary transformation, 

andQ 



H = e^H'e ^ , where H' is the Hamiltonian of [16 



-X 



p + k 



X 



P 



1 + 



2p- k 



k = \k\ 



The second term of Eq. (jT^) is generated by this unitary transformation. We note that the 
non-logarithmic part of SW^a cannot be determined unambiguously, since it mixes with other 
0{1/ c^) operators through IR divergence. Our definition merely represents one possible scheme. 
Only when we add all the contributions to the spectrum at 0{l/c^), the sum is free from IR 
divergence and can be defined unambiguously. Since at present we do not know the full form 
of the Hamiltonian up to (9(l/c'^),f] this problem cannot be circumvented in any case. 



t The operator representation of X reads 



X ^ - 



0{a%). 



^ The ©(l/c^) = 0{a'g) non-logarithmic term of VqcdI^") is not known yet, which would mix with Wna 
through IR divergence. All other 0{l/c^) operators of the Hamiltonian have been identified in Rcf. p6|. 
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/i [GeV] 


1.0 


2.0 


3.0 


4.0 


5.0 


rub [GeV] 


5.458 


5.131 


5.027 


4.969 


4.930 



Table 1: Depencence of the pole mass mb{mb,ag on the scale for (Mz) = 0.1181 and 

rub = 4.190 GeV. 



We have to specify how we treat the pole mass rUb in the operators p'^/rrib, U, Wa, W^a, 
6U and SWna- We express the pole mass in terms of the MS mass frib in the series 
expansion in a^g\^) up to using Eqs. (2), (7) and (8) of [0.0 After that we can, in 

principle, reexpand these operators in a^g\fi), since the pole mass enters the denominators and 
im- There is, however, no known guiding principle how to organize these expansions, since 
such reexpansions cannot be carried out consistently with the expansions in l/c.0 This is in 
contrast to the reexpansion of E^^^{r) = 2mb + VqcdI'"); which has a guiding principle by the 
cancellation of (9(Aqcd) renormalons, although the reexpansion is indeed inconsistent with the 
1/c-expansion. Hence, we keep the pole mass as a function of TiTb, a^g\^i) and fi and do not 
reexpand the operators {p'^/nib, U, Wa, Wna, SU, SWna)', the values of the pole mass are 
shown in Tab. |l] corresponding to a^g\Mz) = 0.1181 and rrib = 4.190 GeV. We will examine 
how the uncertainty of the pole mass of order Aqcd affects our predictions in Sec. 6. 

We take the input for the strong coupling constant as a^g\Mz) and calculate a^g\^2{r)) 
for ii^imp(r) and a^s\fi) for the other operators. We evolve the coupling by solving the 3-loop 
renormalizat ion-group equation numerically and match it to the 4- and 3-flavour couplings 
successively through the matching condition ^^.f\ (Although the 4-loop running of the MS 
coupling constant is available, we consider the 3-loop running more consistent in our analysis, 
which incorporates corrections up to the 2-loop finite part of Vqcd 

Our choice of the zeroth-order Hamiltonian and the way we organize the perturbative expan- 
sion is largely motivated by the success of phenomenological potential models. In fact, the above 
organization of the perturbative expansion follows, to a large extent, the approaches of phe- 
nomenological model analyses, if we identify our H^'"""^ with the non-relativistic Hamiltonian, 
p^/m + Vphcnol?"), used in those analyses. Since i?imp(r) agrees well with typical phenomenolog- 
ical potentials up to an additive constant, we expect that we can make close comparisons with 
phenomenological model analyses, and that eventually the bottomonium spectrum may be re- 
produced with a good accuracy. Note that our zeroth-order quarkonium wave functions, which 
are determined from Hq^'^\ include some of the higher-order corrections in the usual order 
counting of the 1/c-expansion, since Ei^p{r) includes the 0(l/c) and 0(l/c^) corrections to 

^ These formulas were derived originally in ]l|, 0. 

^ For instance, if we organize the perturbative expansions of the energy levels appropriately, renormalons 
contained in them are of order Aq(-,j5, whereas if one expands the levels in 1/c, renormalons become order 
^QCD ^^'^ worsen the perturbative convergence p2| . This may be regarded as an explicit example of problems 
due to an inconsistency between 1 /c power counting and renormalon cancellations (beyond those controlled by 
e-expansion p^). 

II The matching scales are taken as rrib and nic, respectively. 
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/i 


IS 


IP, 


2S 


2P, 


3S 


1.0 


9.476 


9.877 


9.986 


10.186 


10.247 


2.0 


9.498 


9.896 


10.007 


10.203 


10.262 


3.0 


9.505 


9.902 


10.013 


10.209 


10.268 


4.0 


9.509 


9.906 


10.017 


10.212 


10.271 


5.0 


9.512 


9.908 


10.020 


10.214 


10.272 



Table 2: Zeroth-order energies E^^'. The units are GeV. Unless otherwise stated, all tables and 

figures use mf, = 4.190 GeV, a^p{Mz) = 0.1181 and the pole mass mi,{rni„ a^^\fx) , fi) as explained in 
the text. 



the static QCD potential. By the same token, the zeroth-order energy levels are different from 
the Coulomb energy levels, and in particular they depend on the orbital angular momentum /. 

3 Zeroth-order Energy Levels, Wave Functions and Scales 

We show numerical solutions to the zeroth-order Schrodinger equation (|^). Here and hereafter 
we take the input value for the strong coupling constant to be the present world average value 
af\Mz) = 0.1181 p5|. Through Sees. we use the bottom quark MS mass rrib = 4.190 GeV 
taken from [Q. We show the zeroth-order energy levels E^p^ in Tab. ^. The squared radial 
wave functions multiplied by the phase space factor are shown in Fig. ^. Both the energy levels 
and the wave functions depend on the scale /i (only) through the pole mass in p'^/nib in the 
Schrodinger equation. The energy levels Eij^^ vary by about 10 MeV (20 MeV) when /i is varied 
from 2 to 5 GeV (1 to 2 GeV). Nonetheless, if we take the difference of any of the two energy 
levels, the fi dependences cancel mostly. The /i-dependences of the wave functions are fairly 
weak. 

In 1^, Q the scale fi = fi^ for each quarkonium state | ip ) was fixed by minimizing the /z- 
dependence of each energy level calculated in a fixed-order perturbative expansion (minimal- 
sensitivity prescription). The scale fixed in this way turns out to represent the physical size of 
the corresponding quarkonium state fairly well. This was shown in [0 by comparing the scale 
/i^ and the support function defined by 



f^{q)=e{m-q)- drr^ \R^{r)\^ (19) 



where R^{r) is the radial part of the wave function for the state | ^z^). The support function 
represents the support in the momentum-space integral in the calculation of (the major part 
of) the energy level. The Coulomb wave function, evaluated with a^g\n^), was used to compute 
the support function f\p{q); the charm mass effects in loops were not taken into account. Here 
we compare the support function computed with our zeroth-order wave function and the scale 
fixed in both of which include the charm mass effects. The resulting support functions 
for the S states are shown in Fig. |^ and for the P states in Fig. We see that with respect 
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1 2 3 401234560 2 4 6 8 



r [GeV-^] r [GeV"^] r [GeV^^] 




012345 2 4 6 8 



r [GeV-^] r [GeV"^] 

Figure 4: r^|i?^(r)p for different values of // = 1,2,3,4 and 5 GeV, where R^{r) denotes the radial 
part of the zeroth-order energy eigenfunction. The area below each curve is normalized to unity. The 
one curve that visibly differs from the others corresponds to ^ = 1 GeV. Note that the scales differ 
between the plots. 



1 fJ'SS fJ'2S 




0.5 1 1.5 2 2.5 3 

q [GeV] 



Figure 5: Support functions for the S states. The solid curves show the support functions as defined 
in Eq. ([T9|); for comparison of the relevant scales, a^J^\^) is also plotted (dashed curve). Since the 
analysis that we advocate in this work does not attribute scales to the individual states, the scales 
indicated by the dotted lines are taken from Q, Table II. 
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to the treatment in 0, the support functions are shifted shghtly towards higher momentum. 
This is reasonable, since the wave functions calculated here are peaked closer to the origin in 
coordinate space than the Coulomb wave functions (see e.g. Fig. |^ below). Since the scales /i^ 
are located within the IR "tails" of the corresponding support functions f^p{q), we confirm that 
the above interpretaion of the minimal-sensitivity scale /i^ is valid also when we include the 
charm mass effects and the 0(l/c) and 0(l/c^) corrections of Vqcd('^) to the wave functions. 



4 Fine Splittings 

In this section we examine the fine splittings in the bottomonium spectrum within our frame- 
work and compare them with the experimental data as well as with previous theoretical pre- 
dictions in the literature. Within perturbative QCD, it is expected that in principle the fine 
splittings can be computed much more accurately than the individual energy levels. This is 
because the potentials which contain order Aq^^j^r^ renormalons (the static QCD potential and 
Wisia) do not induce the fine splittings, so that these renormalons cancel in the computation of 
the fine splittings. 

The fine splittings (or fine structure) are the level splittings among the states with the 
same principal quantum number n, orbital-angular momentum / > and spin s but with 
different total angular momentum j. Experimentally, the fine splittings have been observed 
among the l^Pj [{n,l,s) = (2,1,1)] and 2^Pj [{n,l,s) = (3,1,1)] states, we therefore examine 
the predictions for these splittings. 

In perturbative QCD the fine splittings are induced first at (9(a;|mfe) = (9(l/c^) by the 
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operator U: 



A4=> . (.e I U I ,e ) ^ ^ x > 0). ,20) 

where 

= -[^ -3^:^^- (2/-l)(2/ + 3) ' ^^^^ 

^ (L.s) = i + -/(/ + !)- + (22) 

This is the fixed-order formula used in [0, ^. In this formula, the Coulomb wave function 
I V'c ) is used to compute the expectation value. The scale dependence of AE^ is large, since 
it is proportional to [q;5'*(/^)]^: three powers of as come from the Coulomb wave function 
{ipc I X I V'c ) 5 and one power comes from the operator U. In 0, |^ the scale is fixed by the 
minimal-sensitivity prescription; c.f. the previous section. 

In our approach, we calculate the fine splittings from AEu = {tp \ {U+6U) {ip)- As compared 
to the fixed-order formula, some of the higher-order corrections are incorporated through the 
wave function | ip ) and the operator 6U. 

In phenomenological approaches, one computes the fine splittings using the same operator 
U but using the wave functions determined from phenomenological potentials: 

A4p'^^"°) = ( Veno I U I Vpheno ) • (23) 



Some of the higher-order corrections to U constrained by the Gromes relation have also 
been incorporated. As an example we compare our results with those of a phenomenological 
model with a Coulomb-plus-linear potential (Cornell potential) 



i?r'''^^'^ = — + K.ornell(r), (24) 

rrib 

^Cornell (r) = -- + ^ (25) 

with K = 0.52 and a = 2.34 GeV"^ 

Tab. ^ compares the expectation values of the operators Uls and Us with respect to the 
Coulomb wave function \ipc)j our zeroth-order wave function | ip ) and the Cornell wave func- 
tion I V'corneii )• The Operators Uls and Us induce the fine splittings through X^s Uls ~ Ds Us, 
where (Xls^Ds) = (-2, -2), (-1, 1), (1, -1/5) for j = 0,1,2, respectively, and (/, s) = (1,1). 
The fine splittings between adjacent levels are therefore Uls — '^ Us for Pi— Pq and 2 Uls+^ Us/5 
for P2 — Pi. The values of a^g\^) in the operators Uls and Us are taken as 0.36 and 0.468 for 
the IPj states and as 0.36 and 0.726 for the 2Pj states. Also the same values of a^s\fi) are 
used for calculating the Coulomb wave functions \ipc)- The first value (0.36) corresponds to 
the value used in the phenomenological analysis [jl[. The latter values (0.468 and 0.726) are 
those for the IPi and 2Pi states which were determined by the minimal-sensitivity prescription 
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as 


Uls 


<l^c) 
Us 


(V^QCD 1 > 

Uls 


<l V^qcd) 
Us 


(V'CorncU > 

Uls 


< V'Cornell) 

Us 


IP 


0.360 
0.468 


2.08 
5.95 


-0.69 
-1.98 


13.71 
17.82 


-4.57 
-5.94 


16.07 
20.89 


-5.36 
-6.96 


2P 


0.360 
0.726 


0.62 
10.22 


-0.21 
-3.41 


8.19 
16.52 


-2.73 
-5.51 


12.92 
26.06 


-4.31 
-8.69 



Table 3: Expectation values (in MeV) of the operators Uls and Us- The mass used in the operators 
is irib = 5.027 GeV corresponding to fni, = 4.190 GeV and /x = 3 GeV. 




2 4 6 8 10 

r [GeV-^] 

Figure 7: Comparison of QCD (soUd), Cornell (dashed) and Coulomb (dash-dotted) IP-wave func- 
tions. For the Coulomb wave function as = 0.36 has been used. 

in 0. Taking into account the numerical values from the table, Uls and Us give roughly the 
same contribution to the Pi-Po-splitting, while the P2-Pi-splitting is dominated by Uls- 

We see that the expectation values of Uls and Us with respect to our wave function | ip ) 
and the Cornell wave function | ^/'corneii ) are very much larger than the expectation values with 
respect to the Coulomb wave function | V'c ) -Q The reason for this behaviour can be understood 
in the following way: Since the potentials P'imp(r) and Vcorneii('") are steeper (i.e. the attractive 
forces are stronger) than the Coulomb potential in the intermediate-distance region^] the wave 
functions are more centered towards the origin for \ip) and | V'Comeii ) than for | ■j/'c ); see Fig. |^ 
which compares the squared radial wave functions multiplied by the phase space factor for the 
IP states. Therefore the wave functions {r\-ip) and (r |'?/'corneii) are enhanced close to the origin 

* This is not necessarily true for the 2P states with as = 0.726, but this can be regarded as originating from 
another effect, which we explain below. Namely, the value as — 0.726 is unrealistically large. 

^ The cancellation of O(Aqcd) renormalons suggests that this behaviour can be understood naturally in 
terms of the QCD force F{r) — — Vqcd('') — —CFOiF{^lr)/r'^ : F{r) becomes more attractive than the Coulomb 
force as r increases due to the running of the F-scheme coupling constant api^/r) ||]. 
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„ 0.014 




12 3 4 

r [GeV-^] 

Figure 8: Comparison of the integrands for (IP | Uls I l-P); the conventions are the same as in Fig. |^. 
as = 0.36 has been used both in the operator Uls and in the determination of the Coulomb wave 
function. Note that the scale of the horizontal axis is different from Fig. ^. 

as compared to the Coulomb wave function. The enhancement factors turn out to be large 
for the P-wave states.0 Since the expectation values of Uls and Us are determined mostly 
by the short-distance behaviour of the wave functions, they are enhanced by large factors. 
This feature can be verified in Fig. H, where we compare the integrands when the expectation 
values of Uls are expressed as integrals over r; the main contributions come from distances 
r < 2 GeV~^. We should stress that the enhancement of the wave functions originates mainly 
from the behaviour of the potentials in the intermediate-distance region and not from the short- 
distance behaviour. This can be seen from the fact that the difference between {ip \ U(^ls,s) I ) 
and (V'Corneii | U(Ls,s) I V'Corneii ) IS much Smaller than the difference between (i/j \ U(^ls,s) I ^ ) and 
{ipc I U(^LS,s) li'c)'- the essential difference between our potential Ei^p{r) and the Cornell po- 
tential resides in the short-distance region. 

We include also the contribution of the 0{l/c^) operator 6U into our prediction of the 
fine splittings. As far as we know, there is no other operator which contributes to the fine 
splittings at 0{l/c^). Considering that our wave function | ip ) includes all the next-to-leading 
order [(9(l/c)] corrections, our prediction for the fine splittings incorporates all the effects up 
to O^agTrib) [0{l/c) relative to the leading C(a|m;,) splittings]. Inclusion of the operator 
6U reduces the scale dependence of our prediction. In Tab. ^ we compare the fine splittings 
calculated from the matrix elements {iplU and from {iplU + 6U The former depend 

■t By way of example, if we squeeze the function r" such that it takes a same value at half distance, i.e. 
r" — > (2r)", then the enhancement factor becomes 2". Thus, the enhancement factor is larger for a larger n. 
For the P-wave states, n = 4, because 2 powers come from the wave function squared and 2 powers come from 
the phase space. 
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asifi) 


IP, 
U 


-IPo 
U + 6U 


IP, 
U 


-IPi 

u + su 


2Pi 
U 


-2Po 
U + 6U 


2P2 
U 


-2Pi 
U + 6U 


1.0 


0.454 


33.3 


22.3 


26.7 


16.2 


20.7 


12.9 


16.6 


9.0 


2.0 


0.301 


22.7 


24.5 


18.2 


18.8 


13.7 


14.4 


11.0 


10.9 


3.0 


0.253 


19.3 


23.3 


15.4 


18.1 


11.5 


13.6 


9.2 


10.5 


4.0 


0.228 


17.5 


22.3 


14.0 


17.4 


10.4 


13.0 


8.3 


10.1 


5.0 


0.212 


16.3 


21.5 


13.0 


16.8 


9.7 


12.6 


7.7 


9.7 


Experiment 


33 


20 


23 


13 



Table 4: Fine splittings computed from {ip\U [ip) and {iplU + 6U ^ is in GeV, the splittings 
are given in MeV. 



on /i rather strongly and are larger for smaller /i, since U is proportional to ^^''(/i).! We see that 
the scale dependence has decreased considerably by the effect of 6U. The scale dependences 
become minimal at ~ 2 GeV for both IPj and 2Pj states. We may try to reproduce the 
splittings calculated with {ijj\U + 6U \ ijj ) a.t these scales by taking an appropriate choice of 
the scales in the lower-order predictions (the splittings calculated with {ip \ U \ ip )) . Then the 
scales become ~ 2 GeV for both the IPj and 2Pj states. All these scales are larger than the 
scales chosen for the respective states in which are fixed by minimizing the scale dependence 
of the individual energy levels, {fi 1 GeV for the IPj states and fi ~ 0.8 GeV for the 2Pj 
states; see Fig. ^) This feature is consistent with a naive expectation: we would expect that 
the latter scales represent the typical scales of the binding energies, or the inverse of the sizes 
of the boundstates (Fig. whereas the former scales represent those probed by the operators 
U or U + 6U, which are larger because the contributions to the matrix elements come from 
shorter distances (Fig. H). 

To our knowledge, so far there has been no systematic argument on the order of renormalons 
contained in the fine splittings. Naively we expect that the largest renormalon contained in 
the calculation of the fine splittings would be of order Aq^p/m^ ~ 1-10 MeV. This estimate is 
based on the Gromes relation which tells us that a part of the operator U is determined 



from the static QCD potential: 

2mlr dr 



L ■ S (26) 



and the fact that dVQC'D{f)/dr contains an (9(AQQj-,r) renormalon. We may take this as an 
order of magnitude estimate of uncertainties of our present predictions for the fine splittings. 
Compared with this error estimate, our predictions of the fine splittings calculated with {ip\U + 
SU lip) in Tab. ^ are in reasonable agreement with the experimental data. 

We are now able to interpret the reasons why the fine splittings of the energy levels computed 
in turned out to be quite small for the IPj states but not so much for the 2Pj states. The first 
point to note is that large enhancement factors are generated by the fact that the quarkonium 
wave functions are more centered toward the origin if we solve the Schrodinger equation with 

^ The scale dependence of the wave functions | through the pole mass is very weak, c.f. Fig. ^. 
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2S - IPi 






u 




Wa 


Wt^A 


Wn A^^Wn A 


1 n 


1 OQ 4 


11 1 


17.6 


-10.6 


-66.5 


25 6 


9 n 




1 5 

XtJ .u 


15.1 


-10.0 


-28.9 


— 1 9 1 

X C/ . X 


s n 

O.KJ 


1111 


1 9 7 

XZj . 1 


13.7 


-9.8 


-20.4 


—90 3 


4.0 


111.4 


11.5 


12.8 


-9.7 


-16.5 


-19.4 


5.0 


111.5 


10.7 


12.2 


-9.7 


-14.2 


-18.3 








3S - 2Pi 






U 




Wa 


Wf^A 


W^A + ^W^A 


1.0 


60.7 


12.6 


9.2 


-0.5 


-29.8 


18.8 


2.0 


59.4 


8.3 


8.2 


-1.5 


-13.0 


-6.8 


3.0 


59.0 


7.0 


7.5 


-1.9 


-9.2 


-8.2 


4.0 


58.8 


6.3 


7.0 


-2.1 


-7.5 


-8.1 


5.0 


58.7 


5.9 


6.7 


-2.2 


-6.5 


-7.8 



Table 5: S-P splittings. The splittings due to the differences in E)^ \ U, U + 6U, Wa, Wna and 
Wna+SWna are given for the 2S'-lPi and 3S'-2Pi splittings. All values are in MeV (/i in GeV). 



Eimp{r) (which we beheve to be more reahstic) rather than with the Coulomb potential. The 
second point is that the natural scales to be chosen for evaluating the expectation values of the 
operator U are larger than the natural scales for the individual energy levels: we confirmed this 
by incorporating the effects of the higher-order operator 6U, and the results are qualitatively 
consistent with a naive expectation. We find that that the first effect overwhelming the second 
one resulted in the quite small splittings among the IPj states in |^. On the other hand, the 
cancellation of the first and the second effect resulted in reasonable sizes of the fine splittings 
for the 2Pj states in that paper, which, in the light of our present observations, may be regarded 
as rather accidental. 

5 S-P Splittings 

In this section we examine the splittings between the S'-wave and P-wave states. In the 
Coulomb spectrum the S'-wave and the P-wave states with the same principal quantum number 
n are degenerate. In perturbative QCD, the splittings are induced by Vqcd('") starting from 
0{l/c) as well as by U, Wa, Wna at C(l/c^). Among these operators, Vqcd('") (after cancelling 
the order Aqcd renormalon) and Wna contain order Aq^j-jr^ renormalons. Therefore, the order 
^qcd'"^ renormalons do not cancel completely in the perturbative computation of the S-P 
splittings. Namely, the theoretical uncertainties of the S-P splittings are expected to be larger 
than those of the fine splittings. 

In Tab. ^ we show the S-P splittings of our zeroth-order energy levels Elp \ which contain 
the effects of Vqcd('") up to 0{l/c^). The contributions of the operators U, Wa and Wna to 
the S-P splittings are also displayed in the same table. (A convenient formula for evaluating 
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2S - IPi 






3S - 2Pi 






+{U,Wa,Wna) 


+{6U,6Wna) 




+{U,Wa,Wna) 


+{6U,6Wna) 


1.0 


109.4 


54.4 


141.9 


60.7 


42.9 


88.2 


2.0 


110.7 


86.8 


96.8 


59.4 


53.2 


59.3 


3.0 


111.1 


93.6 


94.8 


59.0 


54.9 


56.5 


4.0 


111.4 


96.6 


95.1 


58.8 


55.6 


55.7 


5.0 


111.5 


98.3 


95.7 


58.7 


55.9 


55.4 


Exp. 


130 


100 



Table 6: S-P splittings. This table is similar to Tab. H], but here we add up the matrix elements of 
all the operators contributing to the S-P splitting. For comparison, the first of the three columns 
for each level splitting again gives the splitting due to the difference in In the second column 

additionally U, Wa and Wata have been taken into account and finally in the third column 6U and 
6Wna are added. 



the expectation values of Wa is given in the Appendix.) We expect that the contributions 
of these operators would be smaller than the i?^°''-splittings, since the operators are (9(l/c^), 
whereas E^^^ contains the 0(l/c) effects of Vqcd(^)- One sees that this expectation is satisfied 

in most cases, the only operator giving a contribution comparable in magnitude to E^^^ is Wjsia 
for relatively low scales, ~ 1 GeV, where this contribution becomes particularly large. In 
fact, the scale-dependence of the contribution of W^a is large, because Wna is proportional to 
ag'\fiy, whereas U and Wa are proportional only to ag'\^). We see that the scale-dependence 
is reduced by including the effect of the higher-order correction SWna, although Wna + ^W^a 
is still unstable around /z = 1 GeV. If we choose a relatively large scale /i ~ 3 GeV, the contri- 
butions of the higher-order corrections, 6U and 6Wna, become small and the scale dependences 
of the S-P splittings are relatively small; see Tabs. |^ and |[ 

Furthermore, as shown in Fig. ||, the expectation values {ip \ Wna \ ) (V' I Wna + 
6Wna I ^ ) are dominated by short-distance contributions. Therefore, following the same line 
of argument as in the previous section, we expect the natural scales for these contributions to 
be larger than that for the binding energy. 

If we compare our predictions for the S-P splittings with the predictions of the fixed-order 
perturbative expansion in 0, there are two competing effects, just like what we found in the 
case of the fine splittings: these effects are the difference of the wave functions and the difference 
of the scales fj, in the operator. Consequently our predictions for the S-P splittings turn out 
to be larger for a smaller scale fi, are typically larger than the fixed-order predictions for the 
2S'-1P splittings, and are of similar magnitude to the fixed-order predictions for the 3S'-2P 
splittings. Since, however, the S-P splittings are dominated by the 0{l/c) correction from 
VqcDir) [i.e. the contributions from the 0{l/c^) operators are only subleading], the differences 
from the fixed-order predictions are not as pronounced as in the case of the fine splittings. 

As long as we calculate the S-P splittings directly, as done up to here, we see no indication 
of large theoretical uncertainties. However, when we examine the individual energy levels, some 
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12 3 4 

Figure 9: (Absolute value of the) Integrand of the matrix elements of Wna (solid) and SWna (dashed) 
for the IPj states. The graphs peak at rather small distances, indicating that the natural scales for 
the expectation values are larger than those for the boundstate. 





IS 




IPi 


IP2 


2S 


2Po 


2Pi 


2P2 


3S 


u 


13.2 


-25.7 


-6.4 


9.0 


6.3 


-15.4 


-3.8 


5.4 


3.2 


U+5U 


13.2 


-30.7 


-7.4 


10.6 


6.3 


-17.9 


-4.3 


6.2 


3.2 


Wa 


-22.1 


-18.8 


-18.8 


-18.8 


-28.6 


-18.9 


-18.9 


-18.9 


-20.8 


Wna 


-79.1 


-11.9 


-11.9 


-11.9 


-32.3 


-6.4 


-6.4 


-6.4 


-15.6 


Wna+5Wna 


-130.2 


-31.4 


-31.4 


-31.4 


-51.7 


-16.8 


-16.8 


-16.8 


-25.0 




Table 7: Expectation values of the operators f/, U + 6U, Wa, Wna and Wna + ^Wna at the scale 
^ = 3 GeV. 

indications of fairly large uncertainties show up. Let us now investigate this feature. 

Tab. 1^ shows the expectation values of the operators U, U+6U, Wa, Wna, Wna+^Wna for 
all the states and with /i = 3 GeV. We see that the expectation values of Wna and Wna+SWna 
for the IS* state are much larger than what we would expect for 0(l/c^) corrections. Moreover, 
for all the states, the scale-dependences of the expectation values of Wna + ^Wna are large, 
and are comparable to those of Wna', see Tab. ^. (A large part of these scale-dependences 
cancel in the S-P splittings.) The reason for the large scale-dependences is that the non- 
logarithmic term of SWna is very large. Since at present we do not know the full form of the 
Hamiltonian up to 0{l/c^), we cannot draw a definitive conclusion whether this problem of 
large scale-dependence can be remedied. 

The agreement of the S-P splittings between the theoretical predictions and the experimen- 
tal data has improved somewhat from the fixed-order results 0. It is, however, not as good 
as one would naively expect. Namely the differences of the splittings between our predictions 
(for the scale choice fi ~ 2-5 GeV) and the experimental data are larger than the contributions 
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2S 












U 


Wa 


Wna 


Wna+SWna 


U 


U + 6U 


Wa 


Wna 


Wna+5Wna 


1.0 


11.0 


-36.7 


-104.8 


-73.0 


-11.1 


-6.5 


-26.1 


-38.3 


-98.6 


2.0 


7.4 


-30.6 


-45.6 


-63.0 


-7.6 


-7.7 


-20.5 


-16.8 


-44.0 


3.0 


6.3 


-28.6 


-32.3 


-51.7 


-6.4 


-7.4 


-18.8 


-11.9 


-31.4 


4.0 


5.7 


-27.5 


-26.2 


-45.0 


-5.8 


-7.2 


-17.8 


-9.6 


-25.7 


5.0 


5.3 


-26.9 


-22.6 


-40.6 


-5.4 


-7.0 


-17.2 


-8.3 


-22.3 



Table 8: Expectation values of the operators U, U +SU, Wa, Wna and Wna+SWna for the states 
25" and IPi for various scales fi. For the 25 state we leave out U +5U because 6U = for / = 0. 



of the (9(l/c^) operators. As stated in the introduction, it is important to clarify whether the 
level of disagreement is still consistent within perturbative uncertainties. From the examina- 
tion of the contributions to the individual energy levels, we conjecture that there would still 
be large theoretical uncertainties to the S-P splittings, particularly from the contributions of 
Wna, since at present we have no systematic argument on how much of the 0{Aqqj)) renor- 
malons cancel in the S-P splittings. That is, we consider the cancellation of the large and 
scale-dependent contributions from Wna in the S-P splittings to be accidental, unless we find 
a systematic argument in support of it. In this regard, we consider the predictions for the S-P 
splittings much less reliable than those for the fine splittings. 

There are also other indications that the operator Wna becomes a source of instability of 
theoretical predictions. For instance, ref. addresses Wna to be the source of the large 



uncertainties of the cross section for e e~ tt close to threshold. Ref. [29| shows that IR 
logarithms, which are related to higher-order corrections to Wna and Vqcd, generate very 
large corrections to the bottomonium spectrum at (9(a|mb) and consequently cause a large 
uncertainty to it. Furthermore, refs. |32| resum the IR logs; it was found that Wna is 

reduced due to this resummation in the corresponding Wilson coefficients. Ref. P2| claims that 
the resummations lead to a stable theoretical prediction for the top threshold cross section. 



Ref. , however, reported a result which disagrees with that of ^ . 

In this paper we discard the IR logs altogether, which appear first at 0{agmblogas) in the 
spectrum. This is not because we consider them unimportant, but because we believe that we 
are not yet in a position to treat this problem properly. The present treatment of the IR logs in 



|T^, seems to comprise following unsatisfactory aspects in the calculations of the heavy 



quarkonium spectrum and the top threshold cross section (besides the disagreement). Since 
loga^ is not particularly large, it appears that the non-logarithmic terms cannot be neglected in 
comparison to the IR logs. These non-logarithmic terms are not yet fully known at 0{a^mb) 
and beyond, so that we are unable to incorporate them unambiguously. Furthermore, if a 
resummation of the IR logs stabilizes the theoretical predictions substantially, we would like to 
understand the physical meaning behind it, i.e. why the resummation is important. A hint to 
these questions was suggested in it shows that 0{Aqqy)) renormalons will be suppressed 



if we incorporate the offshellness of the quarks, which acts as an IR cutoff in the temporal 
dimension. This applies to the C(Aqcj-,) renormalons contained in Wna and Vqcd (and their 
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/i 


IS 


IPo 


IPi 




2S 


2Po 


2Pi 


2P2 


3S 


1.0 


9.369 


9.818 


9.840 


9.856 


9.982 


10.182 


10.194 


10.203 


10.281 


2.0 


9.424 


9.894 


9.918 


9.936 


10.014 


10.231 


10.245 


10.256 


10.303 


3.0 


9.460 


9.915 


9.938 


9.956 


10.032 


10.246 


10.260 


10.270 


10.315 


4.0 


9.481 


9.927 


9.949 


9.966 


10.043 


10.254 


10.267 


10.277 


10.321 


5.0 


9.494 


9.934 


9.955 


9.972 


10.050 


10.259 


10.271 


10.281 


10.325 


Exp. 


9.460 


9.860 


9.893 


9.913 


10.023 


10.232 


10.255 


10.268 


10.355 



Table 9: Energy levels for all states including E^p, 11+611, Wa, and Wna + ^Wna- The parameters 
used are ag {Mz) = 0.1181 and = 4.234 GeV, this choice of mi, is explained in the text. 

higher-order corrections), and they are closely related to the IR logs and their resummation. 
We suspect that further investigations of all these problems may be a way to clarify and solve 
the problem of the operator W^a we face here. 

6 The Spectrum 

In this section we compare the whole bottomonium spectrum as determined experimentally 
with various theoretical predictions. We list the energy levels numerically in Tab. ^ and show 
the spectrum in Fig. |10[ The levels were calculated according to the framework explained in 
Sec. 1^, using the input parameter af\Mz) = 0.1181 and including the effects of SU and 6Wna- 
We employed the scale choices fi = 1, 2, 3, 4 and 5 GeV. 

In this section we use nib = 4.234 GeV instead of = 4.190 GeV to make the the prediction 
for the IS* state coincide with the experimental value for our favoured scale of fi = 3 GeV.|^ 
We consider this different choice of rrih a more natural way to achieve coincidence with the 
experimental value than a simple shift of the whole spectrum; the numerical difference between 
these two prescriptions is actually very minor. 

In Fig. |TU| we compare our results to the spectrum obtained from the experimental data 
("Exp"), the predictions of the fixed-order perturbative expansions ("BSV")| and the result 
of [0] as a typical prediction of recent phenomenological models ("Eichten-Quigg"). 

One can verify the conclusions of our analysis in the previous two sections: The fine splittings 
and the S-P splittings are larger than those of the fixed-order results for the n = 2 states, 
whereas they are of similar magnitude to the fixed-order results for the n = 3 states. The scale- 
dependence of our predictions originates mostly from the scale-dependence of the operator 
Wna + SWna] c.f. the discussion in the last section. Only the gross level spacings between 

^ The change of TrTi, by this adjustment does not alter the qualitative features of the predictions for the fine 
splittings and the S~P splittings discussed in previous sections. 

II We follow the scheme and the scale-fixing condition A of Sec. 4.2 of except that we use the numerical 
solution to the renormalization-group equation for the strong coupling constant. The results are obtained for 
the input parameters ag {Mz) = 0.1181 and 0.1161; the prediction with the latter input agrees better with the 
experimental data. 
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^ = 1 ^ = 2 /i = 3 /i = 4 /i = 5 Exp Eichten- BSV BSV 
^ ^ ^ ^ Quigg (1161) (1181) 

Figure 10: Comparison of the energy levels obtained with different formalisms. The columns labelled 
/i = 1 through n = 5 show our results, where now a^g\Mz) = 0.1181 and ruh = 4.234 GeV have been 
used to make the IS state coincide with experiment for /i = 3 GeV. The columns for = 1, 2, 4 
and 5 GeV have been shifted to achieve this coincidence. "Exp" shows the experimental values and 
"Eichten-Quigg" those obtained in Q. Finally, "BSV" corresponds to the formalism of the two 
columns represent choices of ai^\Mz) = 0.1161 and 0.1181, respectively. 
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adjacent n's are affected visibly by changes of /i between 2-5 GeV, whereas also the S-P 
splittings vary visibly for a smaller /z between 1-2 GeV. The size of this variation is very large 
considering that the scale dependence is formally an 0{l/c^) = 0(a|mb) effect. The level 
spacings between consecutive n's as well as the S-P splittings increase for smaller /x. We 
regard the large scale-dependence generated by Wna and its higher-order corrections as the 
largest theoretical uncertainty of our prediction. 

Let us note the effects of the operator W^a or W^a + SWna in particular: the level spac- 
ings between consecutive n's are increased, while the S-P splittings are reduced. The reason 
can be understood as follows. The operator Wna{+SWna) generates an attractive potential 
proportional to (with a logarithmic correction) which is particularly strong at short dis- 
tances. Hence, those states which have larger wave functions close to the origin acquire larger 
binding energies. Since the states with lower n have larger amplitude close to the origin and 
therefore acquire larger binding energies, the level spacings between the adjacent n's become 
wider. Since the S states acquire larger binding energies than the P states, the S-P splittings 
becomes narrower. 

Generally, our prediction of the spectrum has a better agreement with the experimental 
data than the fixed-order results. The agreement seems to be better for a larger scale choice, 
which appears reasonable, since the natural scale of the operator Wna + ^W^a^a would be large. 
(The scale dependences due to other effects are much smaller.) The agreement seems to be 
optimal for a scale choice /i ~ 3-4 GeV. We also examined our predictions for different values 
of the input a^g\Mz) within the present world-average values 0.1181 ± 0.0020 ||2^. We find 
that generally all the level spacings and splittings become larger for larger a^g\Mz), since the 
binding energy increases. The widening of the level spacings, however, can be compensated 
largely by choosing a larger value for /i. The spectrum of the phenomenological model still has 
a better agreement with the experimental data, since it includes more parameters which can 
be adjusted. 

In general one should carefully take into account theoretical uncertainties when comparing 
the whole spectrum with the experimental data. Based on the renormalon argument, each 
energy level has (at least) an uncertainty of order Aq^pr^, where r should be taken as a typical 
size of the quarkonium state. Theoretical uncertainties contained in i?imp(r) in the region 
r < riR can be represented typically by these renormalon estimates [|, |[. Ref. estimated 
the uncertainties to be ±(5-30) MeV for the 15* state, ±(20-130) MeV for the n = 2 states, and 
±(40-220) MeV for the n = 3 states. On the other hand, the level spacings (splittings) have 
smaller theoretical uncertainties since these theoretical uncertainties cancel, at least partly, as 
we discussed in the previous two sections. 



Let us estimate the errors of our prediction from other sources. In the left part of Fig. 11 



we show the effects of a variation of the 6-quark pole mass by ±300 MeV from the value listed 
in Tab. |1|. (For the error estimates we set the scale to /i = 3 GeV.) The states with principal 
quantum numbers n = 1 and 2 are shifted up (down) by about 20 MeV when the pole mass 
is shifted down (up) by 300 MeV; for the states with principal quantum number n = 3 this 
variation is about ±15 MeV. Consequently if we compensate the overall shift such that the 15* 
level agrees with the experimental value, only the n = 3 levels vary by about ±5 MeV. 

In the right part of Fig. Ill] we also show the effect of a variation of the long-distance part 
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m6-300 MeV QCD m6+300 MeV iCi QCD 2Ci 



Figure 11: Analysis of various uncertainties: In the left part of the diagram we show the effect of 
a shift of ±300 MeV of the pole mass while keeping the MS mass that enters the potential 
^imp(?') constant. To make the 15 states concidc, wc have shifted the spectrum for m;, — 300 MeV 
(mb+300 MeV) down (up) by about 20 MeV. The right part shows the effect of changing the slope of 
the IR part of i?imp(r') by a factor of 1/2 (2). (The definition of Ci is given in Sec. 2.2.) The effect of 
this change on the IS state is negligible, so that the spectra did not have to be shifted. 
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Figure 12: For comparison we show the bottomonium spectrum with the Wna+^Wna term artificially 



set to and (M^) = 0.1201 and mi) = 4.151 GeV, otherwise the conventions are as in Fig. IC. The 



scale dependence is decreased and the agreement with the experimental spectrum is strongly improved 
with respect to Fig. 



10, indicating a phenomenological preference for a suppression mechanism for 



Wna- 



of the potential E\^^(r) in the zeroth-order Schrodinger equation. We vary the tangent of the 
linear potential at r > rjR by factors of 2 and 1/2. (The first derivative of the potential then 
becomes discontinuous at r = rjR.) The variations of the energy levels are of order ±(2-5) MeV 
for the 2P states and ±15 MeV for the 35" state, while they are smaller than 0.1 MeV (and 
therefore not visible in Fig. |rT]) for the lower states. This can be easily understood because 
only the states with principal quantum number n = 3 have a wave function that extends to 
large enough distances to probe the potential in this region. 

For the sake of comparison, we show in Fig. |1^ the predictions for the bottomonium spectrum 
when we set the Wi^a + ^^na to artificially (in this figure we use ag\Mz) = 0.1201 and 
nib = 4.151 GeV). We see that the scale-dependence reduces (as expected) and there is a much 
better agreement with the experimental data. Phenomenologically this may be taken as an 
indication that the contributions of Wna should be suppressed by some mechanism. 
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7 Conclusions 



We have examined the bottomonium spectrum within a specific framework based on pertur- 
bative QCD. The computation of the individual energy levels includes all the effects up to 
0{asmi)) and that of the fine splittings contains all the effects up to 0{asmb). We have also 
included some important higher-order corrections to the quarkonium wave functions and the 
energy levels through our use of Ei^p{r), which is the characterizing feature of our analysis. 
We have eliminated the scale-dependence from Ei^p{r) by fixing a specific scheme and we have 
discussed scale-depedences originating from the other terms of the Hamiltonian. The agree- 
ment of the fine splittings among the IPj states between the theoretical prediction and the 
experimental data improved drastically as compared to the fixed-order prediction of 0, ^ . We 
find that the centering of the wave functions towards the origin as compared to the Coulomb 
wave functions, due to the strong attractive force in the intermediate-distance region, strongly 
enhances the fine splittings. We also find that, in accord with a naive expectation, the natural 
scales of the fine splittings are larger than those of the boundstates themselves; the latter were 
used in the analyses of [0, ^. The predictions for the fine splittings are stable against the 
variation of the scale /i and are in reasonable agreement with the experimental data both for 
the IPj and 2Pj states. We also examined the S-P level splittings. The agreement with the 
experimental data has improved as compared to the fixed-order results, but the predictions 
are still somewhat smaller than the experimental values. The predictions for these splittings 
are stable against the variation of fi between 2-5 GeV but become unstable for lower scales 
between 1-2 GeV. Natural scales of the S-P splittings are also found to be larger than those of 
the boundstates. On the other hand, the predictions of the level spacings between the adjacent 
n's depend rather strongly on the scale fi. This stems from the large scale dependence of the 
operator Wna and must be regarded as a major source of uncertainties in our predictions. We 
are motivated to choose a relatively large value for the scale fi in view of the dominance of 
short- distance contributions to {ip \ W^a li')- (Other effects are much less scale dependent.) 
If we choose ~ 3-4 GeV, and for a^^\Mz) = 0.1181 and nTf, = 4.234 GeV, the agreement 
between our prediction and the experimental data for the whole bottomonium spectrum is 
fairly good, and is considerably better than the agreement between the fixed-order prediction 
and the experimental data. There seem to be some indications, however, that the contribu- 
tion of the operator W^a reduces the stablity of the theoretical prediction and at the same 
time worsens the agreement between the prediction and the experimental data. At the present 
state we consider our predictions to be consistent with the experimental data within theoretical 
uncertainties. 

Appendix 

We derive a formula which is convenient for evaluating the expectation value of Wa [Eq. (|^)] 
with respect to the eigenstate | ) of Hq^"^^ defined in Eq. (^). We substitute the following 
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operator identities to Wa'- 



'^-^-4n6^'\r), (27) 
^r^^-i?.n.p(r)l. (28) 



— r^r^pip^ = p 
p ^ = nib 



Then one finds 
Wa ^ 



1 






, Cpas 


' 2ml 





Eimpir) 



mf 2mb I r 



(29) 



Hence, the expectation value can be written as 

OjIWaI^) = 



Am, 



- i?i„,p(r) 



+^'(' + i)(^) + ^w3)r^- (30) 

2m^ \r'^/ 

All quantities on the right-hand-side can be evaluated from the radial wave function and the 
energy eigenvalue, which are obtained by solving the Schrodinger equation numerically. 
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Note Added in Proof: 



After this paper was submitted, we received [0] which confirmed in vNRQCD framework the 



result of pNRQCD |Q, and the disagreement we mentioned at the end of Sec. 5 has been 
resolved. 
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